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Abstract 

Inspired by a formal resemblance of certain g-expansions of modular forms and the master 
field formalism of matrix models in terms of Cuntz operators, we construct a Hermitian one- 
matrix model, which we dub the "modular matrix model." Together with an AA = 1 gauge 
theory and a special Calabi-Yau geometry, we find a modular matrix model that naturally 
encodes the Klein elliptic j-invariant, and hence, by Moonshine, the irreducible representations 
of the Fischer-Griess Monster group. 
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4 Prospectus 



Wenn nur ein Traum das Leben ist, 
Warum denn Miih' und Flag'? 
Ich trinke, bis ich nicht mehr kann, 
Den ganzen, lieben Tag! 



1 Introduction 

The recent resurrection of the old matrix model in string theory has invited the community to 
re-examine many intricate properties of random matrix integrals under a new light. The technology 
of using fat matrix Feynman diagrams to cover the Riemann surfaces associated to the string 
worldsheet is once more pertinent. Non-perturbative information about a wide class of A/" = 1 
gauge theories is ascertained from partition functions of bosonic matrix models |21EII3- Moreover, 
the elegant relation between the spectral curve associated to the eigenvalues of a random matrix 
and the special Calabi-Yau geometry that is used to engineer the supersymmetric gauge theory has 
been pointed out. 

A somewhat parallel phenomenon occurred in mathematics some two decades ago. McKay and 
Thompson observed [5] that the coefficients in the Fourier expansion of Klein's modular invariant 
function j{q) are simple linear combinations of the dimensions of the irreducible representations of 
the then freshly-conjectured Monster sporadic group. This observation prompted the "Monstrous 
Moonshine" conjectures [HI that compelled mathematicians to re-investigate the century-old j- 
function with a modern eye. The final tour de force proof of the Moonshine conjectures by Borcherds 
[7j would employ methods that Frenkel, Lepowsky, and Meurman jH] borrowed from vertex algebra 
techniques arising from two-dimensional conformal field theory. 

The j-function is not new to string theory. The modular invariant torus partition function for 
bosonic closed string theory on the 24-dimensional torus obtained by quotienting M^^ on the Leech 
lattice^ is exactly determined by the j-function ||TU] : 

z^^'^ii) = % ir^ = - m + const., m 

^The Leech lattice ALeech is the unique even, self-dual lattice in 24 dimensions with no points of length-squared 
two. See Ref. ^ ^oi details. 
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oo 

where r]{q) = g^/^^ 11 (1 ~ 5'") is the Dedekind 77-function. The central charge of the conformal 

n=l 

field theory (CFT) fixes the constant term in Z{q). The Monster sporadic group is itself the 
automorphism group of (the Z2 orbifold of) this CFT. 

We suggest that the parallels between recent trends in string theory and Moonshine run deeper 
still. Inspired by McKay's often daring and insightful speculations, in this paper we formulate 
an observation that relates the matrix model description of field theories to the Klein invariant 
j-function. We believe that this interrelation is non-trivial and itself hints at intricate and profound 
connexions among modular forms, simple groups. A/" = 1 gauge theories, and Calabi-Yau geometry, 
as well as to random matrix theory. String theory is a unifying principle that blends these disparate 
notions together. 

The observation begins with the form of the g-expansion of the j-function, which is the following: 

j{q) = ^ + bo + biq + b2q^ + .... (2) 

This is reminiscent of an object in matrix models. There is a famous formulation of matrix models 
known as the "Master Field" ^T] , in which all correlation functions in the model can be computed 
without recourse to complicated integrals over infinite dimensional random matrices. For the (Her- 
mitian) one-matrix model, correlators are encoded in the master field M, which is an operator 
acting on a certain Fock space of free probability. This master field is expandable in a basis of 
so-called Cuntz operators {a, a"!"} as 

M = a + tuq + rriia^ + m2(a^)^ + . . . . (3) 

And so an immediate task is evident. What is the matrix model whose master field is the Klein 
invariant j-function? What is its potential, its free-energy, and its density of eigenvalues? Utilizing 
further the Dijkgraaf-Vafa dictionary, what then are the A/" = 1 gauge theory whose Seiberg-Witten 
curve and the geometrically engineerable Calabi-Yau threefold whose special geometry encode this 
modular invariant? All these various players on the stage would naturally encode the g-coefficients 
of j{q), and hence by Moonshine, the irreducible representations of the Monster. 

Of course, the formal similarity we have noted applies if we take the Taylor series for a generic 
function /(g) and compare /(g)/g to the Cuntz-expansion for the Master field. However, the 
resolvent thus defined may not have the correct branch-cut structure to grant us a consistent 
matrix model, even if it could be analytically determined. Happily, the j-function does satisfy these 
requisites. We have chosen to concentrate on j(g) because it is the primitive modular function in 
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that any other meromorphic function that is modular invariant can be constructed from it. To 
consider how modular invariant functions translate to matrix models, it therefore makes sense to 
treat j{q) as the exemplar. Furthermore, the intimate relation between the j-function and the 
Monster group may suggest a new arena for examining the properties of the largest of sporadic 
groups using string theoretic and random matrix techniques. A realization of the Monster in terms 
of geometry and physics is by itself interesting. 

The organization of this article is as follows. Section |21 will briefly review the pieces of our puzzle, 
viz, the elliptic j-function, the master field formulation of the Hermitian one-matrix model, and the 
Dijkgraaf-Vafa correspondence. Next, in Sectional we construct explicitly the matrix model whose 
master field is j{q)- We shall dub this the modular matrix model. The potential for the model will 
enjoy the property that 

V'{z)r.r\q) + -±-^, q = e'-, (4) 

for 2; G [1, 00) C M. We compute the relevant observables in this matrix model. Then, we construct 
the J\f = 1 gauge theory, its full non-perturbative superpotential, and the Calabi-Yau on which the 
theory may be geometrically engineered. In due course, we shall find a hyperelliptic curve which 
encodes the g-coefficients. This is perhaps related to a conjecture of Lian and Yau regarding the 
mirror map of a class of K3 surfaces as l/j{q). We conclude with a discussion of various prospects 
that await our investigation. 

2 Dramatis Personae 

In this section, we introduce three characters who shall play important roles in our narrative, 
viz, the Klein j-invariant, the master field formalism of Hermitian one-matrix models, and the 
Dijkgraaf-Vafa correspondence between supersymmetric gauge theories and matrix models. As all 
these notions figure prominently elsewhere in the literature, this section is included as a concise 
review. 

2.1 The Klein Invariant j-function 

We begin by introducing the most important function invariant under the modular group: the 

Klein invariant function j{q), with q := e^'^*^. The j-function^ is important because it is the unique 

^Though we are following the literature in calling j{q) the Klein invariant j-function, the attribution is in fact 
less straightforward. The earliest reference to the modular j-function of which we are aware is C. Hermite, "Sur 
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modular function in that all meromorphic functions invariant under SL{2,'Z) are meromorphic 
functions of j{q)- In terms of the Jacobi i^-f unctions, the j-function is defined as 



:= 1728 J(v/^), (5) 

The function J{q) is referred to as Klein's absolute invariant. 

The j-function is a modular invariant in the sense that under the SL{2, Z) transformation 

z z' = -, ad — be = 1, a,b,c,dEZ, (8) 

cz + d 

j(q = e^'^*^) = j(q' = e^'^*^'). Thus, to study the behavior of the j-function, it suffices to look at the 
fundamental domain defined by H/SL{2; Z), where H is the upper half 2;-plane, {3m{z) > 0}. The 
j-function maps this domain to the complex numbers: 

jie^""") :n/SL{2;Z) ^C. (9) 

We briefiy relate a number of other properties regarding the j-function jT2] : 

• The j-function is meromorphic on T-C/SL{2; Z). 

• j{q) is an algebraic number, which is rational or integer at special values of q. 

• The g-expansion of the modular invariant j-function is 

j(q) = + 744 + 196884 g + 21493760 g2 + 864299970 (10) 
20245856256 + 333202640600 q^ + 4252023300096 q^ + . . . , 

which is convergent on the fundamental domain. Indeed, on the upper-half plane, 3m.{z) > 0, 
and so g = e27ri<He(z)g-27ram(z) -j^g^^ exponential decay which overwhelms the growth of these 
integer g-coefficients. 



la resolution de I'equation du cinquieme degre," Comptes Rendus, 46 (1858). Dedekind and Kronecker anticipated 
Klein as well. 
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The coefficients in the g-expansion ()10|) are remarkable. They are simple linear combinations 
of the dimensions of the irreducible representations of the Fischer-Griess Monster sporadic 
group (also, the "Friendly Giant"): 



j-function Monster 

196884 = 1 + 196883, 

21493760 = 1 + 196883 + 21296876, 

864299970 = 2 ■ 1 + 2 ■ 196883 + 21296876 + 842609326, 

This beautiful observation due to McKay and Thompson [S], subsequently dubbed by Conway 
and Norton jB] as "Monstrous Moonshine," was proved by Borcherds [Zj. It is truly arresting 
that the most important modular function is related to the largest simple sporadic group. 
(See Ref. [13j for a fascinating tangled history of the Monster sporadic group and the modular 
j-function.) 



2.2 The One-Matrix Model 

Our next ingredient comes from the physics of (0 + l)-dimensional matrix quantum mechanics. We 
briefly recall some relevant facts about the zero- dimensional, one-matrix model. Our discussion 
follows the conventions of the review Ref. ^1]. Let $ be an x iV Hermitian matrix,^ with its 
potential given by V{^). The matrix model is then defined by the partition function Z, which is 
the following random matrix integral: 

- ^o^jp''^('r'-^{-jEn^>))- (11) 

In eq. (jllj) . we have normalized by the volume Vo\{U{N)) of the space of Hermitian matrices, and 
in the second line, we have written the integration over $ in the basis of eigenvalues. The factor 

i<j 

■^Such a Hermitian matrix has N real eigenvalues. One can equally well think of matrices <& in GL{N,€.) with 
eigenvalues distributed along contours in the complex plane rather than along cuts on the real axis. We shall leave 
this for future work |15|. 
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is called the Vandermonde determinant, which in the large- iV limit, induces a repulsion between the 
eigenvalues of The integral in eq. is performed using saddle point methods, obtained from 
the variation of the integrand with respect to a single eigenvalue Aj. This dictates that to 0{1/N), 
that is to say, in the planar limit, 

fEiT^^m). (13) 

The value of the partition function at (9(1/A^) is therefore simply the integrand evaluated at the 
solutions to eq. ()13|) . and hence the planar free energy is 



^ - ]^ log ^ = ]^ (^2 log A(A.) -^Y. ^(^^)) ' (14) 

with Aj the solution of eq. (fTT^ . To solve eq. (fT!?jl . we introduce the following trace, dubbed the 
resolvent R{z) of our matrix model: 



R(z) -^Tr 



z - ^ 

Multiplying both sides of eq. (fT^ by l/{z — Aj) and summing over z, we arrive at the loop equation 
R{zf - 1 R'{z) - i V\z)R(z) + 4- ""'^'[-""^^^^^ = 0. (16) 

The term involving R'{z) drops out to leading order in N. Additionally, in the large- iV limit, the 

density of eigenvalues 

i 

becomes continuous. The resolvent and the normalization conditions then read 

Rlz)=[dX-^, [dXp{X) = l. (18) 



In terms of the eigenvalue density, the saddle point equation (fT!?|l becomes 

2-fdT{^ = -V'{X), (19) 
J A - r g 

where by we mean principal value integration. Likewise the loop equation (jl6p becomes 

R{zr--V\z)R{z)-^^f{z) = 0, (20) 
9 ^9^ 



with 

f{z) := Ag jd\^ {V'{z) - V'{\)) = Ag {v'{z)R{z) - j d\ ^^^-X^ ) ' ^^^^ 

Eq. ()20|1 is now an algebraic equation in the resolvent, known as the spectral curve. 

The key facts that we wish to extract from this discussion are that knowing the functions p{z) or 
R{z) completely determines the partition function Z, and hence the free energy F and subsequently 
all physical observables of the matrix model. Finding these two functions amounts to solving eqs. 
fll8|l and (fTn|l . These we recognize to be Fredholm integral equations of the first kind and Cauchy 
type, and the reader is referred to Refs. [TBl lTTj for a solution thereof (c/. Ref. [18j in this context). 
The idea is that R{z) is always a multi-valued function, for which we describe various branch cuts 
on the real line, the number of which is determined by the number of critical points of V'{z). We 
shall shortly see that the one-cut model will be our chief interest. The discontinuity across the 
branch cut constitutes the solutions of eqs. (fTH|) and (fTUI). In particular, the density of eigenvalues 
captures the multi- valued piece of R{z), 

p{z) = -i, hm {R{z + le) - R{z - ze)) , (22) 

while the potential captures the holomorphic piece, 

- ^V\z) = lim {R{z + ie) + R{z - ie)) . (23) 

2.3 The Master Field Formulation 

At large- A^, there is an especially convenient formulation which computes the physical observables 
of the matrix model. This is the so-called Master Field formalism [llj. The observables are the 
correlators, i.e., the vacuum expectation values (vevs), of the operators Ok = These are 



{Ok) = lim 1 [ Tr O.exp f--Tr V{<^)] 
N^oo N J \ g J 



(24) 



with Z given by eq. (jllll . The information captured in these correlation functions may be extracted 
by computing the trace of an auxiliary matrix known as the master field. We shall quote the relevant 
results and refer the reader to Ref. [12] for further details. 

The master field offers the following remarkable properties that follow from the application of 
free probability theory pUj : 

• There is a free random variable M that provides a realization of the master field. 
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• No integration needs be performed to evaluate the correlators; it suffices to perform algebraic 
manipulations with 

{Ok) = iT[M^] := lim ^(Ti M"). (25) 

• Any correlator in the free probability theory can be reduced to evaluating expectation values 
of powers of M. 

To be more explicit, the master field M defines a linear functional M acting on a Hilbert space. 
In analogy to the familiar ladder operators for harmonic oscillators, a basis for the Hilbert space is 
given by 

(at)"|0). (26) 

The creation operator is the adjoint of a, which annihilates the vacuum |0). The operators a and 
satisfy the Cuntz algebra: 



aa 



t = 1, a^a =|0)(0|, with a|0) = 0. (27) 



This is nothing but the g-deformed Heisenberg algebra, aa'^ — q a^a = 1, for the special case g = 0. 
The operator M associated to the master field M is a function of a and In particular, Voiculescu 
has shown that it is of the special form 

oo 

M{a,a^)=a + ^mn{a^Y (28) 

n=0 

with the {m„} some scalar coefficients of expansion. 

Using eqs. ^I7\ and (PS|). the vevs may be computed by commutations. These will turn out to 
be polynomials in {m„}. The ffist few such Voiculescu polynomials are 

(Ci) = tr[M] := (0|M(a,af)|0) = mo, 

(C2) = tr[M2] := (0|M(a,a^)2|0) =m2 + mi, (29) 
(C3) = tr[M^] := (0|M(a,a^)^|0) = mjj + 3momi + m2. 

At each stage, a new coefficient emerges. Thus, one can recursively compute the coefficients {rrin} 
from the traces {tr[M"']} and hence all physical observables in the matrix model. 

It is convenient to map M(a, a^), which is an operator on the Hilbert space, to a holomorphic 
function in the complex plane: 

^ 00 

ir(z) = - + Vm„;z". (30) 



Z 

n=0 
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Then 

(0|F'(M(a,at))|0) = dzF{K{z)), (31) 

where F is any meromorphic function and C is a contour around the origin. A proof of this 
statement is given in Ref. (TH]. Now, suppose we have found an operator M such that 

tr[M*^] = (0|M(a,a^)^|0), (32) 

for all k. Then the resolvent R{z), which we recall from eq. (jl5p . is simply the generating functional 
of the moments of the matrix M {i.e., the vevs), becomes, due to eq. (jHT|) . 

R^,) = /rfA^ = tr[^^j =f;;.-'=-V[M'=] = f;^-'=~i(0|M(a,at)'=|0) 

Putting K{C) = T, C = K-^ir) =: H{t), we have 

Riz) = --^ / dTH'{T)\og[z-T] = ^ I dT^^ = H{z)=K-\z). (34) 
zTTi Jq 2m J (J T — z 

The point d'appui of the above discussion is that eq. ()34p states that the resolvent is the inverse, 
with respect to composition, of the holomorphic function K{z) obtained from taking the operator 
M{a, a'^) and making the formal replacements a 1/z and (a^)" i-^ z"-. Stated briefly, the resolvent 
is the inverse of the master field. 

Bearing this in mind, we now have an efficient method of generating the Voiculescu polynomials. 
Take a function with a Laurent expansion as follows (which is the form of K{z)), 

f{z) = - + bo + hz + b2Z^ + hz"^ + b^z^ + .... (35) 

z 

Determining the inverse of this function order by order, ^ we find that 



1 ^ &o ^ bo' + bi ^ bo' + 3bobi + b2 ^ &o^ + 6V&i + 2&i^ + 46o&2 + &3 

z z'^ z^ z'^ z^ 

+ ^ + .... (36) 



*This is done as follows. To first order, let / ^(z) = ^ + then z — f ^{f(z)) implies 1 + xz + b^z 



of reads Sfeg + 36i = 2(feQ + 6i) + x, giving x = feg + bi. We may continue ad infinitum. 



1 + 2boz + 0{z'^), giving x = 5o. Iterating to next order, we let / ^(2) = 7 + ^ + ^^^d matching the coefficients 
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We see the Voiculescu polynomials emerging naturally as the coefficients of expansion of the inverse. 
This is how one conveniently generates the expressions in eq. ()29|) . 

Let us entice the reader with a brief preview. We shall explore a particular one-matrix model 
using the technology of the master field. Our interest in this model is closely related to the formal 
substitutions we have made in defining K{z). We observe that the form of the Klein invariant 
j -function is precisely that of the master field of a one-matrix model, and the coefficients of the j- 
function determine the Voiculescu polynomials. We need only make the formal replacement ^ a 
and g" ^— (a^)" in the Laurent expansion of the j-function.^ This paper explores features of the 
matrix model thus obtained. 



2.4 J\f = 1 Field Theory, Matrix Models, and Special Geometry 

Our final ingredient comes from string theory. Recently, Dijkgraaf and Vafa [3 El IH 1221 12n| have 
made the remarkable observation that the holomorphic data, including the superpotential, of an 
= 1 theory in four dimensions can be computed from an auxiliary matrix model. The matrix 
model encodes both perturbative and non-perturbative aspects of the field theory. We shall summa- 
rize a few salient features of the correspondence. For simplicity, we examine an U{n) gauge theory 
with an adjoint chiral multiplet $ and tree-level single-trace superpotential 

p+i . 

iytree($) = J]T^7fcTr $^ (37) 
fc=l 

• An effective low-energy action for $ can be written in terms of the glueball superfield 

5 = Tr (38) 

where Wa is the gauge field strength of the U{n) vector superfield. The theory is treated 
perturbatively in terms of Feynman diagrams, which are considered using 't Hooft's double 
line notation. Only the planar graphs contribute to the free energy! The full non-perturbative 
superpotential in terms of the glueball is 

M^eff(5) = n^FoiS) + S{n\og{S/k^) - 2iiit), (39) 

where Fq{S) is the planar (genus zero) free energy and r = ^ + ^ is the complexified 
coupling. 

'"'Since q^^ ■ q i-^ aa^ = 1 while q ■ q^^ a^a = |0)(0|, this is certainly not a homomorphism. 
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The field theory results are reproduced by an associated matrix model. We promote the chiral 
field $ to an X Hermitian matrix.^ The potential of the matrix model, to be inserted into 
the partition function (|TT|l . is formally the tree- level superpotential ()37|1 . The planar limit in 
the field theory is the same as the large- limit of the matrix model. The glueball superfield 
S is identified with the 't Hooft coupling gN in the matrix model. 

In the original papers [3 El the correspondence between four-dimensional A/" = 1 field theory 
and matrix models is mediated through Calabi-Yau geometry. The theory ()37|) is geometrically 
engineered [211 123 on the (local) Calabi-Yau three-manifold 

{u' + v' + y' + Wl^X^f = /,_i(x)} C C^ (40) 

where /p_i(a;), an order {p — 1) polynomial, is a complex deformation parameter, to be iden- 
tified with its namesake parameter f{z) from eq. (j23) that appears in the spectral curve ()20|). 
Let the volume form on the Calabi-Yau be Q. Then the following special geometry relations, 
as period integrals along compact A-cycles and non-compact 5-cycles, manifestly identify the 
glueball S: 

s,= f n, ur.= ^= [ n. (41) 



The effective superpotential (c/. eq. (j39| ) |26j is 

W,^{S)= f G3Ar] = VA^.n^ + aV^,, Nr.= ! Gs, a:= f Gg, (42) 

where the Calabi-Yau has a three-form fiux = F^^^^ — tH^^^\ with r the type IIB 
axiodilaton. The constant a, independent of the choice of -B-cycle, is identified with the bare 
coupling of the low-energy Yang-Mills theory. 

Now, the non-trivial part of the Calabi-Yau threefold (jiUj) is the hyperelliptic curve 

y' = wLei^y + f,^,ix), (43) 

which encapsulates the Seiberg-Witten geometry [27j and encodes holomorphic information 
about the gauge theory. But this curve is precisely the spectral curve of the corresponding 
matrix model! Geometry thus interpolates between a field theory and its associated matrix 
model. 



^Though we shall sometimes elide this distinction, the rank of the gauge group U{n) and the rank of $ in the 
matrix model are logically separate and should not be confused. In particular, the large- iV limit is not a U{oo) limit 
in the field theory. 
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3 Modular Matrix Models 



Having set the stage for our analyses and the grounds for our speculation, in this section, we will 
introduce our protagonist: the modular matrix models. We have already mentioned in passing at the 
end of ^2.3l the formal resemblance between the Cuntz-expansion of the master field in a Hermitian 
one-matrix model and the g-expansion of the j-function, or indeed of an abundance of modular 
forms. It is now our purpose to perform the analysis in detail and construct the matrix model 
whose master field is the j-function. We shall name this model, and any similar models for other 
modular forms with analogous g-expansion, the "modular matrix model." 

3.1 Seeking the Master 

We recall from eq. (jlUj) that the j-function affords the g-expansion 

I °° 1 
= - + "^n*?" = - + mo + rriiq + ... (44) 

n=0 

with coefficients 

{mo, mi, . . . ,m6, . . .} = 

{744, 196884, 21493760, 864299970, 20245856256, 333202640600, 4252023300096, . . .} . 

(45) 

Next, we recall that the C-number version of the master field, defined as K{z) in eq. (jHUj) is 

oo 

K{z) = - + Yl '^nZ"'. Caveat emptur! We are here establishing an identification not with K{z) as 

n=0 

was done in eq. (I^Uj). but rather with K{q), where q = e^'^*^, i.e., 

^ oo 

j(g)^i^(g) = - + Vm„g". (46) 

This is a simple change of variables. We know that j(g) has the correct expansion that resembles 
the master field and hence can be used in that its coefficients yield the desired vevs. The resolvent, 
according to eq. ()34j) . is the inverse, being careful now to use q = e^'^*^. Therefore, we must 
determine the inverse function of j(e^'^*^) as a function of z. We shall see below that this is a 
well-known function. We summarize by stating that we seek a matrix model whose resolvent is 
given by 

R[z) = r^{e^^''). (47) 
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3.2 The Inverse j-function 

As a brief prelude, we remind the reader that the inverse function to j{q) can readily be determined 
order by order: 

1, , 1 744 750420 872769632 1102652742882 1470561136292880 

J {q) = - + — + — — + i — + g + ^ + •••• 48 

q q'' q-^ q* q'> q^ 

That is, 

r'(^)=E^«^-' ^ = ^"^^ (49) 

fc=0 ^ 

where = afc({"^j}i=o) is the A;-th Voiculescu polynomial in the coefficients of the j-function as 
given in eq. ()48|) and a_i = 1 by definition.'' 

This expansion, to which we shall return, will of course give us no information about the cut- 
structure. Indeed, the previous subsection dictates that it is j~^(e^'^*^) that is required. Now, the 
inverse function of Klein's absolute invariant J(e^'^*^) is well-known; see for example, Ref. (29j. We 
have that 

5 \^ ^ f I 1 1 



s[z) 



2(V3 - 2) r (li)' ,_F, ^; 1 1 _ , (52) 



where 2Fi{a, b; c; z) is the standard hypergeometric function and the range is 

|^"^(^)| > 1, ~ < ^t{J'\z)) < foTzeC. (53) 
Subsequently, we can determine the inverse of the j-function from eq. (0), which we recall to be 

j(e2-*^) = 1728 J(e"*"), (54) 

as 

Henceforth, we send z/1728 — > z and use this rescaled z without ambiguity. 

^The Voiculescu polynomials and thus the coefficients in the series expansion of the inverse of the j'-function are 
related to the Catalan numbers, which are the number of two-dimensional Dyck paths, and to the proper characters 
of the Monster group. We thank J. McKay for pointing out Ref. From the matrix model perspective, the 

Catalan numbers arise from the counting of planar graphs. 
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3.2.1 The Branch Cuts and Multi-Valuedness 



Let us analyze the branch cuts of eq. 2F'i{a,b] c; z) is a single- valued function on C with a 

single branch cut (1, oo), where it is continuous from below, i.e., for z > 1, 



lim 2-Fi {a, b;c;z — ie) = 2-^1 (a, 6; c; z) 
lim 2F1 (a, b;c;z + ie) = 



r(c - a)T{c - b)T{a + b-c+l) 



iFi{a,b;a + b-c+l-l-z) + e^^^^^+'-^^s^ila, b; 



c; z 



(56) 



Therefore, the branch cut we take for 2Fi{a, b;c;l — z), which appears in eq. (|55|) will be (—00, 0). 
Also, we take the branch cut (1, 00) for \/z'^^. In conjunction, therefore, we have that 



r{z) 



and that 





ie 


—ie 


z e (-00,0) 


—e^r{z) + t{z) 


r{z) 


(0,1) 


r{z) 


r{z) 


z G (1, 00) 


r{z) 


r{z) 



, , _ 2^e-f 7r3/2 fl 12 

^^'^ '- r(|) ^^^1^12' 12' 3'^, 



(57) 



s(z 





ie 


—ie 


z G (-00,0) 


—e^s{z) + u{z) 


siz) 


^G (0,1) 


siz) 


s{z) 


z G (1, 00) 


s{z) 


-s{z) 



u{z) :-- 



, r . , 2ze-f 7r3/2 / 7 7 2 

2(V3 - 2)V2 - 1 TTTT — 2Fx — , — ; -; z 

vv ;v 21 ^^2' 12' 3' 



(58) 



Recalling finally the relation ()47|1 . we find that 

• e"3"(s-r)+(t-M) 



i?(;2 + ie) ± i?(2 - k) = < 



-eir(s+r)+(t+«) 
r+s r— s' 



± i — 



z G (— C)0, 0); 

2:G (0,1); 
z G (1, 00). 



(59) 
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3.3 Obtaining the Matrix Model 



Armed with the above mathematical machinery, we are ready to develop the modular matrix model 
associated with the j-function. First, let us find the eigenvalue distribution. We recall from eq. 
fl22|) that the density of eigenvalues is given by the difference of the resolvent along the branch cuts. 
Therefore, we find, upon simplifying eq. that 



A subtlety needs to be pointed out. The eigenvalue density eq. (jUUI) is not real, which is a property 
enjoyed by Hermitian matrix models. We shall apply an easy cure instead. For simplicity, we will 
consider only real densities, i.e., Hermitian matrix models. Although p is complex from (— oo,0), 
it is real from (1, oo). Let us only consider this latter region. It is actually quite remarkable, that 
there is a region such that a Hermitian matrix model could be retrieved. This truncation, i.e., a 
restriction to the space of large- matrices to be integrated, is certainly a common practice, as was 
performed in periodic potentials such as the Gross- Witten model pT] . 

Another caveat is that p{z) is not bounded, which is another desired feature. This could be 
mended by choosing a normalization scheme as follows. We multiply p by a constant A and let p 
be defined from 1 to a such that 



A will depend on a through eq. (jHH). Indeed, we can take a to infinity and A to zero consistently at 
the end of the calculation. We will ascertain the form of A{a) when we impose finitude on physical 
observables later on. In summary then, our density function, together with a plot (recall that we 



^We thank C. Lazaroiu for pointing out certain subtleties involved in using Hermitian matrix models and empha- 
sizing to us the need to work with complex matrices in the context of the Dijkgraaf-Vafa correspondence PU]- We 
could of course study such a model because the techniques that we employ in this paper generalize. We shall leave 
such generalization to future work |15|. 




z e (— oo, 0); 



(0,1); 
z G (1, oo). 



(60) 




(61) 
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have scaled z by 1728), is as follows. 



. 14 



0, 2 < 1; 



P(z) 




z 



(62) 

A benefit of this simplification is that we are now dealing only with a one-cut matrix model instead 
of a two-cut model, thereby reducing the complexity of the computations significantly. 

What then is the original matrix model? We recall from eq. (|23p . that, complementing the 
eigenvalue density, the potential is given by the sum across the branch cuts. Therefore, within our 
chosen region, we have from eq. ()59|) that 

1^,, , ^ ( .T — s .r + s 



-V\z) = Aia) ( t—^ + t'—^ )=A{ r\e'n + -TTl^ > ^ ^ (1, a), (63) 
g \r + sr — sj \ j ^(^e^^^^j j 

which is a purely imaginary function. Consequently, 

Though an analytic computation of the integral seems intractable, V{z) can be numerically deter- 
mined. 

To give a flavor of the form of the potential, we power expand the first few terms as 

V{z) = -gA{a) i (^1.73205 + 0.335531 z^ - 0.0649988 zl - 0.025183 z + 0.0377546 zl 

-0.0151001 zl - 0.00886708 z^ + 0.0159885 - 0.00729369 zl + 0{z^)) . (65) 



The reader may be disturbed by the appearance of the fractional powers, which simply reflect the 
branch point at the origin. Now, since the minimum value of z we take in our model is z = 1, it 
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is more natural to expand about this point instead. Subsequently, we arrive at an ordinary Taylor 
series for the potential which we shall henceforth use: 

/ (z - (z - 1)^ (z - 1)^ 

V{z) = -gA{a)ii2z + 0.0696001 ^ ^ ' - 0.0284334 ^ ^ ^ +0.0168107 ^ ^ ' — 

{z-l)^ iz-V)^ iz-l)'^ 

0.0115703 ^ + 0.00865627 ^ - 0.00682747 ^ + 

5 6 7 

0.00558517 ^ g ' \ +0{z^). (66) 

In possession of the eigenvalue density and the (classical) potential in our d = matrix model, we 
next find the planar free energy, which is the generating function for all the (connected) correlators. 
In terms of the eigenvalue density and the potential, this is computed in the standard way [21] by 
taking the large- iV limit of eq. (fT^ : 



dx p{x)V {x) — j j dx dy p{x) p{y) \og\x — y\ 
j dxp{x) {]^V{x) -logx^ 



Using the above expressions and eq. (|6ip. we can power expand eq. (|67p as a function of the cut-off 
parameter a as: 

F = -0. 108828*^^2 (a- 1)5 - (0.0335904 A + 0.0544111^^^2) (a- l)t + 

(0.0166887A + 0.00414015 z^A^) (a- l)i - (0.00980934^ + 0.00108499*^^2) (a- l)t + 
(0.00643459^ + 0.000424347*^^2) (a- 1)^ - (0.00454071^ + 0.000204799^^^2) (a- 1)^ + 
(0.00337413^ + 0.00011274^^^2) (a - 1)^ + C(a^). (68) 

Now, this is an observable and needs to be finite. From the expansion above we can approximate 
the behavior of our normalization A with respect to the cut-off scale a. We see that 

00 00 
\F\ < J2 + ^^)(« - 1)"^' ^ A^(a - 1)" + C(A2), (69) 

n=3 n=l 

where we have used the fact that all the numerical expansion coefficients Oj in eq. (j68p have modulus 
less than unity. Therefore, we see that the free energy will be bounded if we choose the normalization 

{a-l)-^/y\og{a). (70) 
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3.3.1 Some Salient Features of the Matrix Model 



In the above we have worked in the 2;-variable, in order to conveniently power expand necessary 
physical quantities. We must not forget however, that the interesting number theoretic properties 
of the j-function, and indeed of any modular form, are encoded in the g-expansion. If we take the 
potential of our matrix model, which we recall from eq. (jUHj) is 

n^) = -9A{j-\e^n + -T^). (71) 

and expand about q = exp(27riz), we would find, according to the order by order treatment of eq. 
(gHl), that 

d^^,, //I 744 750420 872769632 1102652742882 1470561136292880 
-i-Viq) = -gA[[- + — + — + + + + . . . 

( 196884 167975456 180592706130 217940004309744 

+ U - 744 + . . 

\ q q'^ q* 

, / 196883 167974712 180591955710 217939131540112 

= -gA g - 744 - 



q q^ q^ q* 



(72) 



Using dz = a g-expansion for the potential of our Hermitian matrix model may subsequently 

be developed: 

../X -9^ f ... 196883 83987356 60197568710 54485001077436 

y{<l) = ^ g- 7441ogg+ + ^ + + J + . 

\ q q^ q-^ q^ 



(73) 

From eq. (j^Uj) . we recall that the expectation values of the matrix model are Voiculescu poly- 
nomials of the expansion coefficients of the master field. The master field we have chosen in our 
theory is the j-function; therefore, the expectation values computed from eq. (j66j) or eq. (j73p will 
be Voiculescu polynomials in the coefficients of the j-function. Indeed, from the coefficients of the 
series expansion of j~^{z) obtained order by order from j{z) in eq. PHj) . we see that 



744 = 744 = mo; 
750420 = 196884 + 744^ = mi + ml; 
872769632 = 21493760 + 3 ■ 744 ■ 196884 + 744^ = ma + 3momi + ml; 

14C 

m3 + 2ml + 4mom2 + Gm^mi + m^; 



20245856256 = 864299970 + 2 ■ 196884^ + 4 ■ 744 ■ 21493760 + 6 ■ 744^ ■ 196884 + 744^ 



(74) 
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which is the expected pattern of Voiculescu polynomials. If we compute the fat Feynman diagrams 
associated with the correlation functions 



using from eq. (j66|l or eq. (f73|l . we shall retrieve precisely the list in eq. (f74|l . 

3.3.2 A Closely Related Matrix Model 

Equally could we have asked ourselves, now that we have a Hermitian one-matrix model whose 
vacuum expectation values are Voiculescu polynomials in the g-coefficients of the j-function, how 
might we establish a similar model whose moments are the g-coefficients themselves? To construct 
this latter model, we need only invert the Voiculescu polynomials. Therefore, we desire a master 
field 



tr[$] = /io = 744; 

tr[$2] = ii^ + iil = 196884; 

tr[$=^] = /i2 + 3/io/ii + /io = 21493760; 

tr[$^] = /is + 2/i? + 4/io;U2 + 6/ig/ii + /i^ = 864299970; 

tr[$^] = /i4 + 5/ii/i2 + 5/io/i3 + 10/io/i? + lO/Xo/^2 + 10/u|j/ii + ^ll = 20245856256; 




(75) 




n=0 



(76) 



such that 



(77) 



Subsequently, the first few coefficients can be iteratively obtained: 



^lQ = 744, /ii = -356652, ^2 = 405710240, /xg = -582814446942, 



(78) 



In other words, we impose that 



= - + Aio + (yUi + fJ'Dq + (yU2 + 3/ioyUi + /io)g^ + . . . , 



(79) 



signifying that 



1 _^ /fO _^ /il + /ip ^ /i2 + 3/io/ii + /ig 

q q^ 



+ .... 



(80) 



21 



Therefore, using eq. (jHUj) . this means that 

^jil/q) = Miqr\ (81) 

and hence we have constructed the inverse for our master field eq. (fTUI) . which, we recall, is the 
resolvent of our desired matrix model. In terms of z, we have that 

M(e2"*")-^ = e-^"*"j(e"^"*')- (82) 

Thus, 

R(z) = Miz)-' = ^log(\j{l/z)] = --^ + ^logj{l/z). (83) 

An analogous route to the preceding sections could be taken. Now, as j{z) has no branch points, 
the only branch cut comes from the logarithm, which extends from to oo. However, j{z) is not 
defined over the real axis, but only on the upper half plane. Therefore, this deceptively simpler 
model escapes the usual analysis in the context of Hermitian matrix models. Our initial choice, 
wherein the vacuum expectation values are the Voiculescu polynomials in the g-coefficients on which 
we have thus far focused, is a more becoming choice. Of course, we could employ the technology of 
more involved {e.g., complex) matrix models on eq. (|83|) . but that is another story, which we shall 
relate another time. This is a lesson in the simple fact that not every meromorphic function can be 
used as a resolvent of a Hermitian matrix model; that the inverse of the j-function could be, and is 
consistently constructable as one, is a pleasant surprise. 

3.3.3 Variations on a Theme by the Master 

The content of Moonshine subsists in the observation that the coefficients in the g-expansion (fTUI) 
encode data about the irreducible representations of the Monster sporadic group. Any resolvent 
that organizes this information faithfully will offer a realization of the Monster group's character 
table in terms of a matrix model and can then be employed to make statements in the language of 
field theory and geometry about the Monster sporadic group itself. While the relations ()47|) and (jH^ 
are particularly natural in the sense that they express the resolvent or the Voiculescu polynomials 
in terms of the j-function in a convenient and direct manner, there are other models that also 
bottle the substance of Moonshine. For example, we could have taken R{z) = exp{2iTij~^{z)) as 
our starting point. That we have made certain choices in our analysis is not to imply that these 
are the only ones that exist. However, the path we have chosen, to consider R{z) = ^^^(e^'^*^), is 
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particularly organic. In this case, we are led to a zero-dimensional, Hermitian one-matrix model 
with a single branch cut that lends itself to simple computations and still distills some essential 
virtues of Moonshine. 



3.4 A Dijkgraaf-Vafa Perspective 
3.4.1 An A/" = 1 Gauge Theory 

Having constructed a Hermitian matrix model which naturally encodes the j'-function, it will be 
expedient to seek yet more physical quantities which bear connexion with j{z). We shall present 
some intriguing speculations in this section. 

As reviewed in §2.4| it is the realization of a recent set of seminal works by Dijkgraaf and Vafa 
[21 ini H] that matrix models of the type described above actually compute the full non-perturbative 
content of an A/" = 1 gauge theory. The tree-level superpotential, according to the correspondence, 
is simply the (classical) potential for the matrix model. We formally replace the Hermitian matrix 
by the single adjoint field $ of the U{n) Af = 1 gauge theory, and obtain the superpotential as 

/ ($-1)2 ($-1)3 ($-1)^ 

Wtreei"^) = ^ A 2$ + 0.0696001 ^ — —^- 0.0284334^ — -^ + 0.0168107^ — — 



($-1)5 ($-1)6 ($-1)^ 
0.0115703 ^ + 0.00865627 ^ - 0.00682747 ' 
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744 (27ri$) + 196883 e~^"** + 83987356 6"*^**+ 



27rz ' 

60197568710 e"^^** + 54485001077436 e"^"^* + ...). (84) 

using eqs. (jUHj) and (f7H|) . The fact that the tree- level potential is non-polynomial need not disturb 
us; the famous Gross- Witten model [21] has a cosine potential, for which a Dijkgraaf-Vafa analysis 
was carried out in Ref. 0. 

We can use the standard Dijkgraaf-Vafa prescription to determine the full non-perturbative 
potential in terms of the glueball condensate 

S = ^Tr W«W^ (85) 

which obtains from the gauge field strength Wa- Accordingly, one identifies S with the 't Hooft 
coupling in the matrix model, and the full superpotential is simply 

W,s^S\ogS + ^Fo (86) 
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using the planar free energy Fq computed in eq. (jbH|) . It is now necessary to restore the powers of N 
in the matrix model, accompanying the powers of the gauge coupling g, on which the cut-parameter 
a depends through normalization. The derivative becomes, from eq. ()fj8j) . 
d 

—Fq = -0.1Q32A2igA'(S?a'(S)(a(S) - 1)^ + 
oS 

{O.OUmii g A{S f a'{S) + A{S) (0.058411 a'(5) + 0.10882 - 0.03359 (a(5) - l)t + 

(0.0023339 i g A{Sf a'{S) + A{S) (0.03539 a'{S) + 0.00217 i A' (S)) - 0.00981 A'{S)) (a(5) - 1)^ + 
+0((a(5)-l)f), (87) 

where a'{S) and A'{S) are derivatives of a and A with respect to S. The functional form of a{S) is 
easily determined by restoring the coupling and further, the appropriate power of A^, into eq. (j61|) . 

viz, 

A.j\,p^,) = l, (88) 
This means that, upon taking the derivative of iS = gN , 

1^/(«) = T 1-T7T^ ' (89) 



dS'^ ' A\ AdS ^ 

where /(a) is the anti-derivative of -^^-z^f- We can then numerically solve this differential 
equation to find the form of a (5) and back-substitute into eq. (jHTj) to determine the effective 
superpotential in terms of the glueball S. The W^^i we obtain at last, has a g-expansion just like 
the potential V did. A pressing question is whether this is a modular form. Unfortunately, due 
to the fact that we cannot integrate eq. (j63|l analytically, the answer to this question eludes us at 
present. As we only have the perturbative expansion of the function, it is difficult to determine 
such properties as modularity. Be that as it may, were W^s to be a modular form, it would provide 
a natural quantum generalization of the j-function. 

3.4.2 A Calabi-Yau Geometry 

Our final venture, having embarked upon a course from the g-expansion of the j-function, shall be 
to the realm of special geometry. We recall from ^2.41 that the Dijkgraaf-Vafa correspondence also 
provides us a singular Calabi-Yau geometry associated to the matrix model which is a generalization 
of a conifold. The singular geometry is that of hyperelliptic curve, given by the spectral curve of 
the eigenvalue densities, transverse to a C^. In particular, embedded in C^, the Calabi-Yau is 

{u^ + t;^ + s{x, y) = 0} C C[x, y, M, v\, (90) 
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with s{x, y) the hypereUiptic spectral curve. For our model, this curve is 

- V'{xf - ^/(x) = 0, (91) 

with f{x) the remainder term given in the spectral curve (j2(jp . It is to be considered as a quantum 
correction (deformation) to the geometry. 

Now our potential, from eq. is V\z) = A (j-\e'^'''') + j^r^^y Therefore, eq. (ED) tells 
us that we have 



z-X 



(92) 



Hence, the spectral curve is fully determined, complete with its quantum correction f{z). 

The above geometrical digression is reminiscent of an earlier proposal circulated amongst math- 
ematicians. It is a known result (c/. Ref. [S2]) that for the family of elliptic curves over given 

by 

y^ = 4x'-^x-^, seF\ (93) 
s — 1 s — 1 

the mirror map z : i— s G at a point of maximal nilpotent monodromy, is precisely 

ziq) = (94) 

the reciprocal of the j-function, namely the inverse of j{q) with respect to multiplication. In fact, 

an old conjecture of Lian and Yau [321 states that for a wide class of K3 surfaces, the g-series of 

the mirror map is a Thompson series Tg{q) for some element g of the Monster sporadic group Wl.^ 

What is the relation between our hypereUiptic curve (PT|) and the family ()93p of elliptic curves? We 

leave this as an open problem. 

Hauptmodul is a function Jg{q) = q^^ + J2'^=i '^n{g)q'\ for g e Although ~ 8 x 10^^ there are "only" 
171 distinct Hauptmoduls of 971. Consider J\[q) — j{q) — 744. The coefficients of the j-function (and also Ji{q)) 
are related to the irreducible representations of DJl. Replacing the m-th coefficient of the Fourier expansion of the 
Hauptmodul Ji (g) by the corresponding representation of 971, we have a formal power series 

Ti{q) = V-iq-^ +0 + Viq + V2q^ + 

where Vr are head representations of 971 {V-i — 1, Vi — 1 ® 196883, etc.) and V — 0^ Vr- Thompson fS| proposed 
that by analogy, for an arbitrary element g G 971, we write the formal power series 

Tg{q) := chv^qig) = chy_, ig)q^^ + + chy, {g)q + chy, {g)q^ + . . . , 

where chy^ {g) is a group character of 971. 
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4 Prospectus 



Grounded in the observation that the formal resemblance between the g-expansion of Klein's ellip- 
tic j-function and the Cuntz-expansion of the master field in a one-matrix model is more than an 
accident, we have in this paper constructed the "modular matrix model." The key property of this 
Hermitian bosonic matrix model is that its master field is by construction j{q). The vacuum ex- 
pectation values, i.e., the observables, in the theory are Voiculescu polynomials in the g-coefficients 
of the j-function. We have computed the eigenvalue density, the resolvent, and the free energy of 
the model. Furthermore, we have established the (classical) potential. Utilizing the Dijkgraaf-Vafa 
correspondence, we have also constructed the associated supersymmetric Af = 1 gauge theory which 
localizes to such a bosonic matrix model. Subsequently, we have computed the full non-perturbative 
superpotential as well as a local Calabi-Yau geometry as a fibration of a hyperelliptic curve on which 
the gauge theory may be geometrically engineered. 

A host of tantalizing questions immediately opens to us. We have much alluded to the complex 
{i.e., GL{N, C)) version of our analysis, which, among others, would give us a matrix model directly 
encoding the g-coefficients of j{q). Does such an extension lead to further intriguing observations? 
Of course, one cannot resist considering other variations of the matrix model. The current resurgence 
of interest in c = 1 matrix models perhaps tempts us to seek the one-dimensional extension to the 
matrix models that we have considered here. Would this quantum mechanics of matrices, by 
promoting $ to in eq. (jHH), offer new perspectives on the j-function? 

It is clear that we should explore the special geometry of the hyperelliptic curve for the modular 
matrix model further. Can one establish a relation between our hyperelliptic spectral curve and 
the Lian-Yau conjecture concerning the family of elliptic K3 surfaces? Now, the latter conjecture 
relates classes of K3 surfaces to the Thompson series of Hauptmoduls. What about our Calabi-Yau 
three-fold, which describes a deformed conifold geometry? How does our geometry, which naturally 
encodes the j-function, relate to these Thompson series? 

Perhaps the most tempting speculation is what we shall call "Quantum Moonshine." Our anal- 
yses have given us a matrix model corresponding to the Klein invariant j-function. Now because of 
Moonshine, our modular matrix model, in encoding the j-function, further encodes the dimensions 
of the irreducible representations of the Monster group. The manifestation of Moonshine, in its 
relation to vertex operators, has already appeared in string theory. In "Beauty and the Beast" 
|lUj . bosonic closed string theory compactified on an orbifold of the Leech lattice with a back- 
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ground Neveu-Schwarz i?-field was shown to have a partition function that encodes the Monster 
representations. How does our story enter into this picture? 

Already in eqs. (|7Hjl and (f7H|) we have encountered some quite large integers that are suggestive 
of this interplay. It was suggested to us^° that the Voiculescu polynomials in the g-coefficients may 
count certain open-string topological invariants associated with the large- geometry. Moreover, 
as we have emphasized, we have performed all our computations in the planar limit, which is to 
say, we have solved the saddle point equation to 0{1/N). The large- iV limit of this matrix model 
allows us, by the Dijkgraaf-Vafa correspondence, to recover details about an A/" = 1 field theory 
whose tree- level superpotential is determined by the resolvent. 

But this is only a leading order result. One can compute the higher 0{1/N^) corrections in 
the matrix model. In the field theory, this corresponds to non-planar diagrams and gravitational 
corrections. Such higher-order terms arise through quantum effects and should proffer a quantum 
modification to the Conway- Norton Moonshine conjectures. In other words, the 0{1/N^) terms, in 
correcting the observables in the matrix model, and hence the Calabi-Yau geometry, should provide a 
natural generalization of the elliptic j-invariant and hence provide a quantum version of the Monster 
group. What interpretation do the quantum corrections have in terms of the Monster sporadic 
group? Is this a quantum deformation to the presentation of the Monster? We have mentioned in 
^3.41 that the full non-perturbative superpotential might be a modular form, which would be yet 
another way of finding quantum corrections to Moonshine. The Dijkgraaf-Vafa program therefore 
would offer not only a perturbative window into non-perturbative physics, but also a perturbative 
window into number theory. 

Finally, we see that as the form of g-expansion is generic to modular forms, what would an 
analogous analysis to ours give for one of a myriad other modular functions? What number theoretic 
and geometric data would these encode? We have perhaps raised as many questions as we have 
presented answers, and this is our hope. It is our desire that "modular matrix models" should shed 
light into various hitherto unexplored corners of string theory, number theory, and geometry. 



'We thank A. Iqbal for pointing this out to us. 
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